MPI-Ph/2001-005 
Revised version, 
February 7, 2008 



Tests of Power Corrections for 
Event Shapes 
in e~^e~ Annihilation 

P.A. Movilla Fernandez, S. Bethke, O. Biebel, S. Kluth 

Max-Planck-Institut fiir Physik, 
Werner-Heisenberg-Institut 
Fohringer Ring 6 
80805 Munich, Germany 



Abstract 

A study of perturbative QCD calculations combined with power corrections to model 
hadronisation effects is presented. The QCD predictions are fitted to differential 
distributions and mean values of event shape observables measured in e"'"e~ annihi- 
lation at centre-of-mass energies from y/s = 14 to 189 GeV. We investigate the event 
shape observables thrust, heavy jet mass, C-parameter, total and wide jet broad- 
ening and differential 2-jet rate and observe a good description of the data by the 
QCD predictions. The strong coupling constant as(-^zo) ^i^d the free parameter of 
the power correction calculations ao(2 GeV) are measured to be 

as(Mzo) = 0.117ll[JS and ao{2 GeY) = O.bUt^ofi ■ 

The predicted universality of is confirmed within the uncertainties of the mea- 
surements. 
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1 Introduction 



The study of hadronic final states in e"'"e~ annihilation allows precise tests of the theory of 
strong interaction, Quantum Chromo Dynamics (QCD), using event shape observables for 
the analysis of hadronic events. For event shape observables perturbative QCD predictions 
in 0(ag) and in some cases also in the next-to- leading- logarithm-approximation (NLLA) 
are available. The various experiments at the PETRA, PEP, TRISTAN, LEP and SLC 
colliders collected a large amount of data at centre-of-mass (cms) energies ^/s = 14 to 
189 GeV which can be used to make precise quantitative tests of QCD. 

Precision tests of perturbative QCD from hadronic event shapes require a solid un- 
derstanding of the transition from the perturbatively accessible partons to the observed 
hadrons, the hadronisation process. Hadronisation effects cannot be described directly by 
perturbative QCD and are usually estimated by phenomenological hadronisation models 
available from Monte Carlo event generators, e.g. JETSET/PYTHIA [||, HERWIG § 
or ARIADNE \§. 

Alternatively, analytical approaches are pursued in order to deduce as much infor- 
mation as possible about hadronisation from the perturbative theory. Hadronisation 
contributions to event shape observables evolve like reciprocal powers of the hard inter- 
action scale a/s (power corrections) An analytic model by Dokshitzer, Marchesini 
and Webber (DMW) of hadronisation valid for some event shape observables derives the 
structure of the power corrections from perturbative QCD. The model assumes that the 
strong coupling remains finite at low energy scales where simple perturbative calculations 
break down [^0]. The model parametrises the magnitude of non-perturbative effects by 
introducing moments of the running strong coupling as as parameters to be determined 
by experiment. 

Several experimental tests of power corrections in the DMW model with differential 
distributions or 1st moments (mean values) of event shape observables measured in e^e~ 
annihilation have been done [[7|-p!6[]. In the present paper we test power corrections in the 
DMW model to the differential distributions and mean values of event shape observables 
measured in e"'"e~ annihilation experiments at -y/s = 14 to 189 GeV. 

We use resummed 0(ag)+NLLA QCD calculations combined with power corrections 
to fit the event shape distributions with as(Mzo) and the non-perturbative parameter ao 
as free parameters. In the case of the mean values 0{ag) calculations together with power 
corrections are fitted to the data. We investigate the prediction of the DMW model that 
the non-perturbative parameter does not depend on the specific event shape observable, 
i.e. that it is universal. 

Section § starts with an overview of the observables and briefly explains the theoretical 
predictions. The data used in our study and the fit results are presented in section 0. In 
section § we give a summary and draw conclusions from our results. 

2 QCD Predictions 
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2.1 Event Shape Observables 

We employ the differential distributions and mean values of the event shape observables 
thrust, heavy jet mass, C-parameter, total and wide jet broadening. The mean value of 
the differential 2-jet rate based on the Durham algorithm is used as well. The definitions 
of these observables are given in the following: 



Thrust T The thrust value is given by the expression [T^, 18 

where pi are the momentum vectors of the particles in an event. The thrust axis Ut 
is the vector n which maximises the expression in parentheses. We use 1 — T in this 
analysis, because in this from the distribution is comparable to those of the other 
observables. A plane perpendicular to nx through the origin divides the event into 
two hemispheres Hi and H2 which are used in the defintion of heavy jet mass and 
the jet broadening observables below. 

Heavy Jet Mass The invariant mass of all particles contained in hemisphere 
Hi or H2 is calculated |[T^. The observable Mh is defined by 

Mh = max(Mi, M2) /^/s . 

Some experiments use the definition = max(Mi, M2)^/s where in 0{as) we 
have the relation 1 — T = M^. 

Jet Broadening The jet broadening measures are calculated by ||20|| : 

Bk 





Pi X UtI 


2Ei \P^ 





for each hemisphere Hk, k = 1,2. The total jet broadening is given Bj^ — Bi -\- B2 
and the wide jet broadening is defined by i?w = max(i?i, B2). 



C-parameter The C-parameter is defined as pT],p^ 

C = 3(AiA2 + A2A3 + A3A1) 
where A^, k = 1,2,3, are the eigenvalues of the momentum tensor 

Ei \Pi\ 

Differential 2-jet rate The differential 2-jet rate is determined using the Durham jet 
finding algorithm [^]. In this algorithm the quantity yij = 2mm{Ef,Ej){l — 
cos 9ij) / E^^g, E^is = J2kEk, is computed for all pairs of (pseudo-) particles with 
energies Ei, Ej in the event. The pair with the smallest i/ij is combined into a 
pseudo particle by adding the 4-vectors and the procedure is repeated until all 
Vij > Vent- The value of ycut where the number of jets in an event changes from 
three to two is called j/3. The differential 2-jet rate is defined by the differential 
distribution of 1/3. 
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2.2 Perturbative QCD predictions 



We use QCD predictions in 0(ag) matched with resummed NLLA calculations in our 
analysis of differential distributions P^-pB|. In the NLLA the cumulative distribution 



R{y) = /o' l/crtot(dcr/d?/')dy' of an observable y is considered. The NLLA is valid in re- 
gions of phase space where y is small, i.e. where the emission of multiple soft gluons 
from a system of approximately back-to-back and hard quarks dominates (2-jet region). 
QCD predictions in 0{al) are expected to be valid in regions of phase space where em- 
mission of a single hard gluon dominates (3-jet region). We choose to combine the C(as) 
with the NLLA calculations with the ln(i?)-matching scheme, because it has theoretical 
advantages and is also preferred in experimental analysis H,P7|-PB[1. Other match- 

ing schemes exist and will be considered in the study of systematic uncertainties, see 
section below. 

The complete perturbative QCD prediction renormalised at the scale /iR for a cumula- 
tive distribution RpT{y) using the ln(i?)-matching scheme takes the following form [p4| , |27|| : 



InRpT^y) = Lgi{Las{fi^)) + g2{Las{^lR)) (1) 
-{GuL + Gi2i:2)as(/iR) - (^22^^ + G23L^)aUfiR) 

+A{y)as{fiK) + i^B{y) - 2A{y) - ^^(y)') &Uf^R) , 



where L = \n{l/y) and ds = as/(2vr). The functions gi and g2 represent the all-orders 
resummations of leading and subleading logarithmic terms, respectively, and the Gnm 
coefficients are given e.g. in The coefficient functions A{y) and B{y) are defined by 
A{y) = jQ{dA/dy')dy' and B{y) = /q (di?/d?/')d?/', respectively. The differential distribu- 
tions dA/dy and dB/dy are obtained by integration of the 0{al) QCD matrix elements 



using the program EVENT2 . The prediction is normalised to the total hadronic cross 
section evaluated in C(as). 

The renormalisation scale /xr is identified with the cms energy ^/s = Q of the measure- 
ment. The dependence of the perturbative QCD predictions on the renormalisation scale 
is studied by introducing the renormalisation scale parameter = /xr/Q and making the 
replacements of [E^, equation (23). 



The mean values of event shape observable distributions are defined by 

ry,iiax 1 do" , 

{y) = / y—-T-^y , (2) 

Jo (Ttot dy 

where ymax is the largest possible value of the observable y (kinematic hmit). The per- 
turbative QCD prediction of mean values {y)prp in C(ag) is given by 



{y)p^ = Ay&siijR) + i^By + (7r/5oln(x^) - l)2Ayj ag(/iR) (3) 

where = (33 — 2rif)/(127r) with the number of active quark fiavours nf = 5 at the 
cms energies considered here. The 0{as) and 0(«g) coefficients Ay and By are taken 
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from [^. Calculations of mean values in NLLA are not yet available, because the NLLA 
predictions diverge for very small values of y and do not vanish at the kinematic limits 
Umax of the observables. 

2.3 Power Corrections 

Non-perturbative effects to event shape observables are calculated in the DMW model 
as contributions from gluon radiation at low energy scales where perturbative evolution 
of the strong coupling breaks down. The location of the divergence in the perturbative 
evolution of as, known as the Landau pole, is given by Arg^ ^ 200 MeV in the MS 
renormalisation scheme. The model assumes that the physical strong coupling remains 
finite at scales around and below the Landau pole. A new free non-perturbative parameter 

1 /"w 

ao(/Wi) = — / as{k)dk (4) 
Hi Jo 

is introduced to parametrise the unknown behaviour of as{Q) below the so-called infrared 
matching scale /ij. The non-perturbative and the perturbative evolution of the strong 
coupling are merged at the scale /ij which is generally taken to be 2 GeV P^ . 

The power corrections are calculated including two loop corrections for the differential 
distributions of the event shape observables considered here |T^,^. The effect of hadro- 



nisation on the distribution obtained from experimental data is described by a shift of 
the perturbative prediction away from the 2-jet region: 

^JL^'l!^(y-VD,) , (5) 

dy dy 

where y = 1 — T, M^, C, and i?w The factor V depends on non-perturbative 



parameter and is predicted to be universal [36|: 



with the colour factor Cp = 4/3. The factor K = (67/18 — tt'^/6)Ca — (5/9)nf originates 
from the choice of the MS renormalisation scheme. The Milan factor A4 accounts for 



two- loop effects and its numerical value is 1.49 ||3^. The theoretical uncertainty of Ai is 
about 20% due to missing higher order corrections The quantity Dy depends on the 



observable 13, 36 



L>i_T = 2, Dm^ = 1, Dc = 3n, (7) 



Db = ablia^ + Ff,{b,as{bQ)) , b = BT,By^, 0^^ = 1,0^^ = ^- 

A simple shift is expected for 1 — T, and C whereas for the jet broadening variables 
Bt and i?w an additional squeeze^ of the distribution is predicted. The more complex 

^The term "squeeze" refers to the form of event shape distributions which are more peaked in pertur- 
bative predictions compared to the predictions including hadronisation effects. 
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behaviour for the jet broadening observables calculated in [13| is related to the inter- 
dependence of non-perturbative and perturbative effects which cannot be neglected for 
these observables. The power corrections for St and i?w predict in addition to the shift 
an increasing squeeze with decreasing cms energies. The necessity of an additional non- 



perturbative squeeze of the jet broadening distributions was already pointed out in [|T5 

The power corrections for mean values of 1 — T, and C are obtained by taking the 
first moment of equation and read: 

{y) = {y)pT: + VDy . (8) 



In the case of mean values of and -Bw the predictions from [131 are used. For the ob- 



servable 1/3 the leading power correction is expected to be of the type l/Q^ or (In Q)/Q^ 
but the corresponding coefficients are not yet calculated. 

3 Analysis of the Data 
3.1 Data Sets 

Experimental data below the TP peak are provided by the experiments of the PETRA 
(12 to 47 GeV, about 50000 events in total), PEP (29 GeV, about 28000 events in total) 
and TRISTAN (55 to 58 GeV, about 1200 events in total) colliders. Data around the 
TP resonance are from the four LEP experiments with (9(10^) events per experiment 
and from SLD (about 40000 events) while data above the T^ are exclusively from the 
LEP experiments with (9(10^) events per experiment from = 133 to 183 GeV and 
C(IO^) per experiment at y/s = 189 GeV. For event shape distributions table |1| gives the 
references and also the ranges considered in the fits (see section ^]2| below). For mean 
values we consider published data available in the energy range of 13 up to 189 GeV 
B, P^O|-p^, |29|, pO|, |3^, |38|-[55[1 . All data used in this study are corrected for the limited 



resolution and acceptance of the detectors and event selection criteria and are published 
with statistical and experimental systematic uncertainties. 

3.2 Fit Procedure 

The standard analyses use the entire data sets as described above. For the perturbative 
predictions of differential distributions we employ the matched resummed (9(ag)-|-NLLA 
QCD prediction given by equation ([l|) while the power corrections are implemented accord- 
ing to equation (||). For mean values the 0{al) perturbative prediction from equation (||) 
is used combined with the power corrections according to equation (||). 

For each observable we perform simultaneous x^"fits with Q;g(Mzo) and ao(2 GeV) as 
free parameters. The strong coupling ag(Mzo) is evolved to the renormalisation scale 
/^R = x^Q with Q = ^/s of a given event shape distribution or mean value using the two- 
loop formula for the running coupling ||56[. The is defined by = J2i{{di ~ ti)/cri)^ 



where di is the value of measurement i, ti is the corresponding theoretical prediction and 
(Ji is the quadratic sum of statistical and experimental systematic uncertainties of di. 
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The fit ranges for fits of event sliape distibutions are defined individually for each 
cms energy such that the 2-jet region of the distribution is exploited as far as possible. 
The fit ranges are limited by the demands i) that the of extreme bins do not 
contribute substantially to the total of the distribution, ii) that the perturbative QCD 
prediction is reliable, and iii) that the power corrections are under control. Requirement 
iii) is checked by monitoring the ratio of the theoretical predictions without and with 
power corrections, respectively, using the fit results for as(Mzo) and ao(2 GeV). Figure ^ 
(solid lines) presents these ratios for 1 — T, Mh, M^, i?T, -Bw and C at y/s = 35 GeV. The 
fit ranges are chosen such that regions of rapidly varying power corrections are excluded. 
The chosen fit ranges are listed in table |I[ 



3.3 Effects of bb Events at low y/s 

The presence of events from the reaction e"'"e~ — bb at low cms energies y/s can distort 
the event shape distributions, because the effects of weak decays of heavy B-hadrons on 
the topology of hadronic events cannot anymore be neglected. An additional problem 
arises from comparing QCD calculations based on massless quarks with data containing 
massive quarks at i/i close to the production threshold. 

At y/s <^ Mgo bb events constitute about 9% of the total event samples. Ideally one 
would correct the data experimentally by identifying bb events and removing them from 
the sample. However, since we have only published event shape data without information 
on specific quark flavours we resort to a correction based on Monte Carlo simulations. 
We generate samples of 10^ events at each ^/s with the JETSET 7.4 program with 



the parameter set given in ^Tj. For each event shape observable we build the ratio 
of distributions calculated with u, d, s and c quark events to those calculated with all 
events. This ratio is multiplied with the bin contents of the data to obtain corrected 
distributions. This procedure is applied to all data at ^/s < M^o. We correct the mean 
values for the contribution from b quarks using exactly the same procedure as for the 
correction of the differential distributions. It was verified that the simulation provides an 
adequate description of the data at all values of ^/s < Mgo. Figure |^ shows the ratio of 
distributions of 1 — T calculated using u, d, s, and c quark events or all events obtained 
at ^/s = 14 to 55 GeV as an example. The correction is reasonable within the fit ranges 
at ^/s > 14 GeV while at ^/s = 14 GeV the correction is a large effect. 

Systematic effects due to uncertainties in the Monte Carlo parameters are expected 
to be small for the ratio except for those parameters which only affect the bb events 
in the samples. The most important such parameter is the value of e-b in the Peterson 



fragmentation function |Q which controls the fragmentation of b quarks in the simulation. 
Threshold effects on the fraction of bb events at low ^/s which depend on the value of the 
b-quark mass in the simulation are found to be negligible for the fit results. 
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3.4 Fit and Systematic Uncertainties 

We consider the following for both fits to differential distributions and to mean values 
unless specified otherwise: 

Fit error The fit errors for as{Mzo) and ao(2 GeV) are taken from the diagonal elements 
of the error matrix after the fit has converged. 

Renormalisation scale Systematic uncertainties from perturbation theory are assessed 
by varying between 0.5 and 2.0. The changes in the fit results w.r.t. the standard 
results are taken as asymmetric systematic uncertainties. In the case that both 
deviations have the same sign the larger one defines a symmetric uncertainty. 

Matching scheme As a further systematic check in the analysis of differential dis- 
tributions we use different matching schemes, namely the modified \n{R)- and 
i?-matching schemes, to combine the 0(ag) with the resummed NLLA calcula- 
tions . A possible matching scheme uncertainty is defined by the larger deviation 



caused by using \n{R)- or /^-matching. 

Power corrections Uncertainties due to the power corrections come from the choice of 
the value of ni and from the theoretical uncertainty of the Milan factor Ai. We 
vary fii by ±1 GeV and M. by ±20% and take in both cases changes of the fit 
results w.r.t. the standard results as asymmetric systematic uncertainties. No error 
contribution from the variation of fii is assigned to ao(/^i)) because setting /ii to a 
different value corresponds to a redefinition of ao(yUi). 

Fragmentation of b quarks The standard analysis is carried out with corrected data at 
^/s < Mzo based on the JETSET tuning of |^ as explained in section The value 



of the JETSET parameter is varied around its central value = 0.0038 ±0.0010 
by adding or subtracting its error and the analysis including correction of the data 
at ^/s < Mzo is repeated. Deviations w.r.t. the standard results are considered as 
asymmetric uncertainties. 

Experimental uncertainties We examine the dependence of the results on the input 
data taken for the fits in several ways: 

1. We perform the fit without the LEP/SLC data at ^/s ~ Mzo. 

2. In the case of fits to distributions the fits are repeated using seperately either 
the data below or above the Z° peak. This also checks for possible higher 
order non-perturbative contributions to the power corrections. Such tests are 
impractical with mean values, because the sensitivity to the power corrections is 
reduced when only restricted ranges in ^/s are used, in particular for y/s > Mzo. 

3. A further source of systematic uncertainty in the analysis of differential distri- 
butions only comes from the choice of the fit ranges. The lower and the upper 
edges of the fit ranges of all distributions of a given observable are varied in 
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both directions by one bin. We take the largest of the four deviations w.r.t. 
the standard result systematic uncertainty. 

For distributions the largest deviation from 1. and 2. w.r.t. the standard results is 
added in quadrature with the uncertainty from 3. and the result defines a symmetric 
systematic uncertainty. For mean values only the deviation from 1. is taken to define 
a symmetric systematic uncertainty. 

The total errors of the standard results are defined as the quadratic sum of the fit 
errors, the renormalisation scale uncertainty, the power correction and the experimental 
uncertainties. In the case of distributions the larger of the matching scheme and the renor- 
malisation scale uncertainties is included in the total error and the fit range uncertainty 
is added as well. 

3.5 Results of Fits to Event Shape Distributions 

Our standard results for an observable are obtained with = 1 and /ij = 2 GeV. The 
results from the fits are listed in tables || and ^ for as(Mzo) and ao, respectively. The 
signed values indicate the direction in which as(Mzo) and ao(2 GeV) changed w.r.t. the 
standard analysis when systematic effects are studied. The fit curves of the standard 
results and the corresponding experimental data for 1 — T, Mh or M^, C, and 5w 
are shown in figures ^ to Values of x^/d.o.f. and of the correlation coefficients from the 
fits are given in table ^. 

We generally observe a good agreement of the predictions with the data within the fit 
ranges, as indicated by the values of x^/d.o.f.The two fit parameters are anticorrelated 
with a correlation coefficient pat — —80%. The results for as(Mzo) and ao{2 GeV) are 
consistent with each other within the total errors in the case of 1 — T, C and B^. 

The agreement between data for _Bw at ^/s < Mzo and the QCD prediction is not 
as good as with the other observables, see also table ^. The value for as (Mzo) obtained 
for i?vv is about 15% smaller than the values from the other observables. In the case of 
i?w the QCD prediction in the 3-jet regions tends to lie above the data leading to smaller 
values of as (Mzo) in the fit. Fitting only data at ^/s < M^o leads to a significant deviation 
of as(Mzo) w.r.t. the standard result, see table |. This may indicate that yi-dependent 
non-perturbative effects are not fully modelled by the calculations for i?w 

In order to disentangle the different contributions to this effect we performed a fit of 
the 0{al) QCD prediction for i?w combined with power corrections with as (Mzo), and 
ao(2 GeV) as free parameters and obtained as(Mzo) = 0.106 ± 0.001, = 0.10 ± 0.02 
and ao(2 GeV) = 0.65 ± 0.03 with x^/d.o.f. = 0.5. Since the value for as(Mzo) is 
comparatively small and the value for ao(2 GeV) is comparatively large we conclude that 
both the 0(ag)+NLLA perturbative predictions and the power correction calculations 
contribute to the small values of as (Mzo) and large values of ao(2 GeV) observed in 
the standard fits. Small values of as (Mzo) in fits with i?w using C(as)+NLLA QCD 
calculations have also been observed in P,P|,|77|, ^|3^ -pT|. 
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The results for ao(2 GeV) are consistent with each other within about two standard 
deviations of the total errors; in particular the values for ao(2 GeV) from Mh or 
and -Bw are approximately 25% larger than the other results. We note the coincidence 
that Mh and i?w are calculated using only the hemispheres containing more invariant 
mass or transverse momentum, respectively. We conclude that ao(2 GeV) is approxiately 
universal within the total uncertainties of the individual measurements. The results for 
ao(2 GeV) are also consistent with earlier measurements 0,|TT],|T^. 

The values of as{Mzo) obtained from the fits are systematically lower than correspond- 
ing results which use the same (9(ag)+NLLA perturbative predictions but apply Monte 
Carlo corrections instead of power corrections P|,P, P7| , |^|3D| -pT|. From the experimental 
point of view there is a lucid explanation for the differences between the as(Mzo) results 
based on the power corrections and those based on Monte Carlo corrections |15| . The lat- 
ter induce a stronger squeeze to all distributions than the power corrections which simply 
predict a shift for 1 — T, and C without any presence of a squeeze. Although the 



situation improved for the jet broadening observables due to the revised calculations [|T3 
the effect of the squeeze remains below the expectation of the Monte Carlo hadronisation 
models. As a consequence the two-parameter fit favours smaller values for as(Mzo) in 
order to make the predicted shape more peaked in the 2-jet region and hence chooses 
large values for ao in order to compensate the shift of the distribution towards the 2-jet 
region. 

Figure |l] compares hadronisation corrections as predicted by power corrections and by 
the JETSET Monte Carlo program as used in section |3.3| . The hadronisation corrections 
from the Monte Carlo simulation are defined as the ratio of distributions calculated using 
the partons left at the end of the parton shower (parton-level) and the stable particles 
(r > 300 ps) after hadronisation and decays (hadron-level) . In all cases and in particular 
for Mh and Sw the Monte Carlo corrections increase the slopes of the perturbative 
predictions more than the power corrections leading to larger values of q;s(Mzo) in fits of 
the predictions to the data. 

It turns out that the power correction uncertainties for as(Mzo) are negligible for each 
observable while there are significant power correction uncertainties for ao- We conclude 
that q;s(Mzo) is mainly constrained by the perturbative prediction rather than by the 
power correction contributions while is mostly determined by the power correction 
calculations. The strong dependence of on M. is due to the anticorrelation seen in 
equation (|^). 

The total errors are generally dominated by the theoretical uncertainties. We observe 
significant variations of as(^zo) from i?T and -Bw when considering only data with ^/s < 
Mzo in the fits; for By^ this variation is the largest contribution to the total error of 
as{Mzo). For the jet broadening observables the matching scheme uncertainty is larger 
than the renormalisation scale uncertainty for ao(2 GeV) and also for ag(Mzo) in the 
case of St- We also notice that the ao(2 GeV) results from 5w and Mh or M^ have the 
largest power correction uncertainties. 
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3.6 Results of Fits to Mean Values 



The main fits to the mean values of 1 — T, M^, B^, i?w and C are performed with 
as(Mzo) and q;o(2 GeV) as free parameters using = 1 and fii = 2 GeV. In figure |^ the 
results of the fits and the corresponding perturbative contribution (|/)prp of equation (||) 
are shown. The size of the power suppressed contribution is the difference between the 
dashed and the sohd curves in figure 0. Tables || and ^ list the results of the fits and the 
variations found from the studies of systematic uncertainties. 

We find that the fitted QCD predictions describe the data well with x^/d.o.f. ~ 1. 
The fit results for as(^zo) and ao(2 GeV) for all observables are consistent with each 
other within their total errors and have correlation coefficients pat — —90%. The results 
are also generally consistent with the results from fits to distributions. We note that in 
contrast to the analysis of distributions the results for as(Mzo) and Q!o(2 GeV) from 
and -Bw are compatible with results from the other observables. 

For the observable (7/3) we investigated power corrections of the form (lnQ)/Q^, 
1/Q, (In Q)/Q and omitting power correction terms, introducing ai{fii) = (l//ii)^ ■ Jq^ k ■ 
as{k)dk as the second and an unknown coefficient Dy.^ as the third fit parameter |]35| , |6^ . 
All fits yielded x^/d-o.f. ^ 1. For the 1/Q and [In Q)/Q corrections large values for 



as(Mzo) were obtained which are incompatible with the world average within the 
fit errors. Corrections of the and (lnQ)/Q^ type gave 1 to 2% increased values of 

as(^zo) and a value of ai(2GeV) = 0.25 ± 0.03(fit). The results for Dy^ were —0.2 and 
—0.4 for l/Q"^ and {\nQ)/Q'^, respectively, but also consistent with zero within the fit 
errors. We conclude that the data prefer one of these latter types of power corrections 
although the size of the correction is too small to be determined from the available data. 
The smallness of the fitted Dy.^ coefficient justifies to neglect any power correction for 
fits of the (j/3) data and we only quote the result of such fits in table |[ The result for 
as(Mzo) from (7/3) is also in good agreement with the world average value of the strong 
coupling. 

3.7 Combination of Individual Results 

The individual results are combined to single values for q;s(Mzo) and ao{2 GeV), respec- 
tively, following the procedure described in [|, 0]. The combination is done separately 
for the results from event shape distributions or mean values. A weighted average of 
the individual results is calculated with the square of the reciprocal total errors used 
as the weights. For each of the systematic checks the weighted averages for as(Mzo) 
and ao(2 GeV) are also determined and the total error of the weighted average is calcu- 



lated exactly as described in section |3]^. This procedure accounts for correlations of the 
systematic errors. 

We obtain as combined results from the analysis of distributions 

as(Mzo) = 0.1111 ± 0.0004(fit) ±0.0020(syst.)+!|;J|!3i(theo.) 

ao(2 GeV) = 0.579 ± 0.005(fit) ± 0.011 (syst.)+[5;j7J(theo.) . 
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The error contributions refer to the fit error (fit), the variations of the input data sets and 
the fit ranges (syst.) and the variations of the matching scheme, renormahsation scale, 
Milan factor and (theo.). The total correlation coefficient is estimated as p = —0.16 
(see below). The small value for as(^zo) compared to the world average as(^zo) = 
0.1181 ± 0.0034 ||6^ is caused by the small values of as(Mzo) from Mr or and in 
particular B^j. If the results from B^^ are omitted from the weighted averages the results 
become as(Mzo) = 0.1126t|]:[]|j| and ao(2 GeV) = 0.558lg:[!^^. This value for as(Mzo) is 



in better agreement with the world average and with other measurements P, p7| , |53| , |59| -pl[] 
while the result for ao(2 GeV) changes only slightly. 

The results from the study of mean values based on (1 — T), (My), (-Bt), (-Bw) and 
(C) are 

as(Mzo) = 0.1187 ± 0.0014(fit) ±0.0001(syst.)+J5;S(theo.) 
ao(2 GeV) = 0.485 ± 0.013(fit) ± O.001(syst.)+|5;|;2(theo.) . 



The error contributions are defined as explained above for distributions. The estimate of 
the total correlation coefficient is p = +0.17. The values for q;s(Mzo) and ao(2 GeV) are 
in reasonable agreement with the results from distributions, especially when the average 
of results from distributions is calculated without the values from i?w 

Figures ^ a) and b) present the results for as(Mzo) and ao(2 GeV) with error ellipses 
based on the total errors. The correlation coefficients are determined as follows. For ev- 
ery systematic test a covariance matrix is constructed using the systematic uncertainties, 
symmetrised if neccessary, and a correlation coefficient psyst- In cases where the correla- 
tion between systematic deviations of as(Mzo) and ao(2 GeV) of a given systematic test 
has the same sign as the correlation p^t from the fit result we set psyst = Pm- In cases 
where the signs from the correlations from the standard fit and the systematic test are 
opposite we set psyst = +1 or —1 taking the sign from the correlation of the systematic 
test. All covariance matrices are added and the result defines the error ellipsis. Tables ^ 
and ^ show the correlation coefficients obtained with this procedure. The correlation coef- 
ficients of the averages are calculated as the weighted averages of the individual correlation 
cofficients using the products of the individual total errors for as(Mzo) and ao(2 GeV) 
as weights. The figure illustrates that the individual results for as{Mzo) and ao(2 GeV) 
from distributions and mean values are compatible with each other and with the averages 
within the total errors. We consider this as a confirmation of the predicted universality 
of the non-perturbative parameter a^. 

Finally we combine the results for as(Mzo) from the analysis of distributions, mean 
values of 1 — T, M^, i?T, -Bw and C and from {y^) by calculating error weighted averages 
based on the symmetrised total errors. As errors of the final combined results the smaller 
of the errors of the individual results are chosen. The final result for as{Mzo) is 

as(Mzo) = 0.117ll0S . 

The final result for q;o(2 GeV) is obtained by combining the results from distributions 
and mean values again quoting the smaller of the total errors of the individual results as 
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the final errors: 

ao(2 GeV) = 0.5131°:°^^ _ 

The total correlation coefficient is again estimated as a weighted average of the correlation 
coefficients from the combined results from distributions and mean values, respectively, 
yielding p = +0.07. 



4 Summary and Conclusions 

The analytic treatment of non-perturbative effects to event shape observables in e^e~ 
annihilation based on power corrections was examined. We tested predictions for the 
differential distributions and mean values of the observables 1 — T, Mh or M^, B^, i?W; 
C and 1/3, respectively. For this test a large amount of event shape data collected by 
several experiments over a range of e"'"e~ annihilation energies from ^/s = 14 to 189 GeV 
was considered. 

Fits of perturbative QCD predictions combined with power corrections to distribu- 
tions and mean values of event shape observables were performed with the strong coupling 
as(Mzo) and the non-perturbative parameter ao(/ii) as free parameters. The good quality 
of the fits with x^/d.o.f. ~ 1 supports the predicted 1/Q evolution of the power correc- 
tions. The results for as{Mzo) and ao{2 GeV) from distributions are more consistent with 
each other compared to previous studies [|1^,|16| due to the improved predictions of power 
corrections to the jet broadening variables. However, we still observe a large deviation of 
the as(Mzo) results obtained from i?w from those extracted from the other observables. 
We conjecture that this discrepancy is a combined effect of the perturbative (9(ag)+NLLA 
predictions and the power correction calculations for this observable. The individual re- 
sults for as(Mzo) from all observables are observed to be systematically smaller than the 
corresponding results in P, p7|, pB|, |5U|-pT[] , which use Monte Carlo hadronisation models. 



This observation may be related to the different amounts of squeeze of the distributions 
predicted by both types of hadronisation model. 

We obtain as combined results for the strong coupling constant and the non-perturbative 
parameter: 

as(Mzo) = 0.1171^Ho2o 
ao(2GeV) = 0.5131°;°^^ _ 



It should be noted that the values for as{Mzo) and ao(2 GeV) from i?w are only compat- 
ible with the combined result and with the values of as{Mzo) from the other observables 
within about two standard deviations of the total errors. 

The average value for ao(2 GeV) is in good agreement with previous results P, |35| , |6^ . 
The scatter of ao{2 GeV) values derived from 1 — T, and C is covered by the expected 
theoretical uncertainty of the Milan factor of about 20% [BB . 



Since this value is representative of the individual results within the errors we consider 
this as a confirmation of the universality of ao as predicted by the DMW model. However, 
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the results from Mh or and 5w from distributions indicate that uncalculated higher 
orders may contribute significantly to the non-perturbative corrections. 
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Tables 





Experiment 


1 -T 


Mh, M2 






C 


1 so 


L3 |11| 
OPAL ^ 


U.Uzo — U.oU 
0.03-0.30 


n nQ n 1 fi 
U.Uo — U.io 

0.14-0.45 


n c\P. n OP. 
U.UO — U.zO 

0.05-0.25 


nn/ic; nine; 
U.U40 — U.iyo 

0.04-0.20 


n 1 n n c; 
U. iU — U.OO 

0.08-0.60 


183 


DEI 
L3 

OVI 


.P 
11 
VL 


HI [12| 

H 


0.03-0.28 
U.Uzo — U.oU 
0.03-0.30 


0.03-0.20 

n nQ n 1 fi 
U.Uo — U.io 

0.14-0.45 


0.05-0.24 

n nfi n 9r 
U.UO — U.ZO 

0.05-0.25 


0.03-0.20 

n n/1 K nine; 
U.U40 — U. iyo 

0.04-0.20 


0.08-0.72 

n 1 n n 7n 
U. iU — U. ( U 

0.08-0.60 


172 


DEI 

L3 [ 
OVI 


.P 
11 
^L 


HI [12| 
11 


0.04-0.32 

0.025 - 0.30 

n riQ n Qn 
U.Uo — U.oU 


0.04-0.20 

0.03-0.18 

n 1 /I r\ AK 
U.i4 — U.40 


0.06-0.21 

0.06-0.26 

n nc; n oc^ 
U.UO — U.zo 


0.04-0.17 

0.045 - 0.195 

n n/i n on 
U.U4 — U.zU 


0.08-0.64 

0.10-0.70 

n nfi n fin 
U.Uo — U.OU 


161 


DEI 
L3 

OVI 


.P 
11 
^L 


HI 

m 


12 




0.04-0.32 
0.05-0.30 
0.03-0.30 


0.04-0.20 
0.03-0.18 
0.14-0.45 


0.06-0.21 
0.06-0.26 
0.05-0.25 


0.04-0.17 
0.045 - 0.195 
0.04-0.20 


0.08-0.64 
0.10-0.70 
0.08-0.60 


133 


ALI 
DEI 
L3 
OPi 


IV 
.V 

VL 


H |3 
HI| 

m 


9| 
12 




0.04-0.30 

n n/1 n qo 
U.U4 — u.oz 

n nc; n oc; 
U.UO — u.zo 

0.03-0.30 


n n/i n on 
U.U4 — u.zu 

n nQ c\ ^K 

U.UO — U.iO 

0.14-0.45 


n nfi n 9i 

U.UO — U.Zi 

n nfi n 9r 
U.UO — u.zo 

0.05-0.25 


n n/i n 1 7 

U.U4 — U. i ( 

n n/1 nine; 
U.U40 — u. iyo 

0.04-0.20 


n nfi n fi/i 

U.Uo — U.04 

n T n n 7n 
U. iU — U. ( u 

0.08-0.60 




ALI 

DEI 

L3 

OPi 

SLL 


IV 
.V 
44 
VL 
) 1 


H |6 
HI 1 

m 


5| 
42 

5/ 


1 


n nfi n Qn 
U.UO — U.oU 

0.06-0.30 

n nfic; n qq 
U.UDO — u.oo 

0.06-0.33 

n nfi n Q9 
u.uu — u.oz 


n nQc; n 1 R 
U.Uoo — U.io 

0.04-0.16 

n nQQ n i fiQ 
u.uoy — u.ioO 

0.20 - 0.40 
n nzL n 1 fi 

U.U4 — u.io 


0.09-0.27 

0.09-0.26 

n nfi n 9fi 
u.uo — u.zo 


0.06-0.17 

0.06-0.18 

n nfi n 9n 
u.uo — u.zu 


n 1 fi n 79 
U.iO — U. ( Z 

0.16-0.72 

n 1 fi n 7n 
U.iO — U. ( u 

0.16-0.64 
n 1 fi n fi4 

U.iO — U.04 




AMY 


Din n '^n 

U. lU — U.OU 










AA 


JADE i,i 
TASSO |5[ 


u.uu — U.OU 

0.06-0.32 


n 99 _ n zL9 

0.06-0.16 


n 1 n _ n 9zL 

U. lU — u.z^ 


n nfi _ n 1 fi 

u.uu — U. lU 


n 1 fi _ n 79 

U. lO — U. 1 z 


35 


JADE ||,[ 
TASSO 5 


1 


0.06-0.30 
0.06-0.32 


0.22-0.38 
0.06-0.16 


0.10 - 0.24 


0.06-0.16 


0.20-0.72 


29 


HRS |43[ 
MARKII 1 


49 




0.10-0.325 
0.10-0.32 










22 


TASSO i; 


1 


0.10-0.32 


0.06-0.18 








14 


TASSO i; 


1 


0.12 - 0.32 


0.10-0.18 









Table 1: The sources of the data and the fit ranges for the observables 1 — T, Mh or M^, 
i?T, -Bw and C are shown. The cms energy ^/s at which the experiments analysed their 
data is given in GeV. The observable Mh is used only by OPAL and JADE while M^ is 
used by the other experiments. 
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1 - T 


Mh,M2 

n5 hi 






C 




0.1173 


0.1105 


0.1114 


0982 


0.1133 


fit error 


itU.UUUO 


itU.UUUD 


itU.UUUD 


itU.UUUO 


itU.UUU4 


mod. ln[K) 


1 n nn 1 o 

+O.UUio 


+U.UUU5 


+U.UU5o 


+u.uuiy 


1 n nn 1 /I 

+U.UUi4 


mod. R 


-0.0010 


-0.0005 


-0.0021 


-0.0012 


-0.0005 


Xfj^ — 0.5 


-0.0041 


-0.0023 


-0.0039 


-0.0011 


-0.0039 


x^, = 2.0 


+0.0055 


+0.0037 


+0.0050 


+0.0023 


+0.0052 


M - 20% 


+0.0003 


+0.0002 


+0.0001 


< 0.0001 


< 0.0001 


M + 20% 


-0.0003 


-0.0002 


-0.0001 


< 0.0001 


+0.0001 


= 1 GcV 


+0.0007 


+0.0004 


+0.0003 


+0.0001 


+0.0001 


yi^i = 3 GeV 


-0.0006 


-0.0004 


-0.0002 


-0.0001 


-0.0001 




±0.0002 


< 0.0001 


< 0.0001 


±0.0001 


< 0.0001 


> Mzo 


+0.0023 


-0.0022 


-0.0005 


-0.0008 


< 0.0001 


< Mzo 


-0.0022 


-0.0006 


-0.0032 


-0.0068 


-0.0004 




-0.0026 


-0.0014 


-0.0006 


-0.0002 


-0.0018 


fit range 


±0.0013 


±0.0010 


±0.0009 


±0.0014 


±0.0005 


total error 


+0.0063 
-0.0051 


+0.0045 
-0.0034 


+0.0063 
-0.0063 


+0.0073 
-0.0072 


+0.0056 
-0.0044 



Table 2: Values of as (Mzo) are shown derived from fits of resummed 0(Q;g)+NLLA 

QCD predictions combined with power corrections to distributions of the event shape 
observables 1 — T, Mu or M^, Bx, -Bw and C. In addition, the statistical and systematic 
uncertainties are given. Signed values indicate the direction in which as (Mzo) changed 
with respect to the standard analysis. 
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1 - T 


Mh,M2 

XT 7 J-£ 






C 


an(2 GeV) 


4Q2 




655 


787 


507 


fit error 


1 n nno 


1 n ni 1 


1 n ni n 


1 n ni 
itU.UlO 


1 n nnt^ 

itU.UUO 


mod. m[K) 


— O.Olo 


—0.013 


—0.070 


— 0.0o4 


—0.049 


lliULL. I V 


-l-fl DDR 


-i-D nni 


-i-D m Q 








—0 012 


—0 021 


—0 010 


—0 038 

\J . WOO 


—0 015 


r =20 


+0.009 


+0.014 


+0.007 


+0.027 


+0.011 


M — 20% 


-1-0 063 


-1-0 1 51 


-1-0 113 


-1-0 1 55 


-1-0 077 


M + 20% 


-0.042 


-0.101 


-0.075 


-0.103 


-0.051 




±0.003 


±0.001 


< 0.001 


±0.001 


< 0.001 


v/i > Mzo 


-0.050 


+0.073 


+0.012 


+0.039 


< 0.001 


< Mzo 


+0.028 


-0.018 


-0.005 


< 0.001 


+0.003 




+0.029 


-0.014 


-0.010 


-0.040 


+0.008 


fit range 


±0.019 


±0.018 


±0.018 


±0.024 


±0.003 


total error 


+0.084 
-0.070 


+0.170 
-0.128 


+0.135 
-0.105 


+0.175 
-0.131 


+0.092 
-0.071 



Table 3: Values of ao are shown derived from fits of resummed (9(ag)+NLLA QCD pre- 
dictions combined with power corrections to distributions of the event shape observables 
1 — T, Mh or M^, i?T, -Bw and C. In addition, the statistical and systematic uncertainties 
are given. Signed values indicate the direction in which ao changed with respect to the 
standard analysis. 





1 -T 


Mh, M2 




Bw 


C 


standard fit 


172/263 


137/161 


91.9/159 


96.1/132 


150/208 


fit correlation 


-0.88 


-0.75 


-0.85 


-0.81 


-0.82 


total correlation 


-0.17 


-0.10 


-0.47 


-0.32 


-0.17 


> Mzo 


73.5/131 


43.3/97 


67.4/115 


50.1/100 


95.2/134 


v^ = Mzo 


43.2/59 


68.0/33 


16.1/28 


22.9/20 


44.2/49 


v/i < Mzo 


55.3/69 


25.4/27 


8.4/12 


23.1/8 


10.6/21 



Table 4: The values of x^/d.o.f. and the correlation coefficients between q;s(Mzo) and 
aQ are shown for the standard fit in the first two rows. The third row shows the total 
correlation coefficients between as (Mzo) and ao including effects of systematic variations 
of the analysis (see section |3.7| for details). The other rows present values of x^/d.o.f. 
obtained from fits with subsets of the data with ^/s > Mzo, ^/s = Mzo or ^/s < Mzo, 
respectively. 
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(1-T> 


(Ml) 


(Bt) 




(C) 


(2/3) 


as(-^zo) 


0.1217 


0.1165 


0.1205 


0.1178 


0.1218 


0.1199 


fit error 


±0.0014 


±0.0016 


±0.0015 


±0.0015 


±0.0014 


±0.0008 


xVd-o.f. 


50.1/41 


24.0/35 


23.7/28 


10.4/29 


18.4/26 


13.2/15 


fit corr. 


-0.89 


-0.89 


-0.91 


-0.94 


-0.93 


n.a. 


total corr. 


0.23 


0.18 


-0.09 


-0.47 


0.22 


n.a. 


= 0.5 


-0.0048 


-0.0026 


-0.0037 


+0.0017 


-0.0045 


-0.0039 


= 2.0 


+0.0059 


+0.0037 


+0.0048 


+0.0003 


+0.0056 


+0.0050 


M - 20% 


+0.0020 


+0.0011 


+0.0014 


+0.0009 


+0.0020 


n.a. 


M + 20% 


—0.0016 


—0.0009 


—0.0012 


—0.0008 


—0.0016 


n.a. 


/ii = 1 GeV 


+0.0009 


+0.0005 


+0.0006 


+0.0004 


+0.0009 


n.a. 


Hi = 3 GeV 


-0.0009 


-0.0005 


-0.0006 


-0.0004 


-0.0008 


n.a. 


Eh ± Icr 


±0.0002 


±0.0001 


±0.0001 


±0.0001 


±0.0002 


< 0.0001 




+0.0008 


-0.0020 


-0.0012 


+0.0005 


+0.0015 


+0.0030 


total error 


+0.0065 
-0.0054 


+0.0047 
-0.0038 


+0.0054 
-0.0044 


+0.0025 
-0.0025 


+0.0064 
-0.0053 


+0.0059 
-0.0050 



Table 5: Values of as (Mzo) are shown from fits of 0(as) QCD predictions combined 
with power corrections to mean values of 1 — T, M^, B^, 5w, C and y^. Statistical and 
systematic uncertainties are also given. Signs indicate the direction in which 0:3(^20) 
changes w.r.t. the standard analysis. The renormalisation and infrared scale uncertainties 
are added asymmetrically to the errors of as (Mzo). 





(1-T) 








(C) 


ao(2 GeV) 


0.528 


0.663 


0.445 


0.425 


0.461 


fit error 


±0.015 


±0.021 


±0.020 


±0.029 


±0.013 


Xn — 0.5 


+0.002 


+0.010 


+0.021 


+0.118 


+0.004 


Xf, = 2.0 


-0.001 


-0.003 


-0.014 


-0.046 


-0.002 


M - 20% 


+0.072 


+0.107 


+0.055 


+0.048 


+0.054 


M + 20% 


-0.049 


-0.072 


-0.037 


-0.032 


-0.037 


ffb ± Ic 


±0.002 


±0.007 


±0.002 


±0.002 


±0.002 


v/i ^ Mzo 


-0.003 


+0.017 


+0.007 


-0.005 


-0.006 


total error 


+0.074 
-0.051 


+0.111 
-0.078 


+0.063 
-0.045 


+0.131 
-0.063 


+0.056 
-0.040 



Table 6: Values of ao are shown from fits of 0{a^) QCD predictions combined with 
power corrections to mean values of 1 — T, M^, B^, B^^ and C . Statistical and systematic 
uncertainties are also given. Signs indicate the direction in which ao changes w.r.t. the 
standard analysis. 
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Figure 1: The figure presents hadronisation correction factors at ^/s — 35 GeV estimated 
using power corrections (solid lines) or using the JETSET Monte Carlo program (dashed 
lines). The hadronisation corrections for power corrections are given by the ratio of the 
perturbative QCD prediction over the same prediction combined with power corrections 
using the fitted values of Q!s(-^zo) and q;o(2 GeV). The Monte Carlo hadronisation cor- 
rections are given by the ratio of distributions calculated at the parton- and hadron-level, 
respectively. 
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Figure 2: The figure presents ratios of distributions of 1 — T calculated using u, d, s, 
and c quarks events or all events using Monte Carlo simulation. The different line types 
indicate the cms energy at which the Monte Carlo simulation was run. 
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Figure 3: Scaled distributions for 1 — T measured at ^/s = 14 to 189 GeV. The error 
bars indicate the total errors of the data points. The solid lines show the result of the 
simultaneous fit of as{Mzo) and ccq using resummed 0(Q;g)+NLLA QCD predictions with 
the ln(R)-matching combined with power corrections. The dotted lines represent an 
extrapolation of the fit result. 
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Figure 4: Scaled distributions for Mh and measured at ^/s = 14 to 189 GeV. The 
error bars indicate the total errors of the data points. The solid lines show the result of 
the simultaneous fit of q;s(^-^z") and ao using resummed ©(aD+NLLA QCD predictions 
with the In (R) -matching combined with power corrections. The dotted lines represent an 
extrapolation of the fit result. 
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Figure 5: Scaled distributions for and measured at ^/s = 35 to 189 GeV. The 
error bars indicate the total errors of the data points. The solid lines show the result of 
the simultaneous fit of Q!s(Mzo) and cto using resummed ©(aD+NLLA QCD predictions 
with the In (R) -matching combined with power corrections. The dotted lines represent an 
extrapolation of the fit result. 
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Figure 6: Scaled distributions for C measured at ^/s = 35 to 189 GeV. The error 
bars indicate the total errors of the data points. The solid lines show the result of the 
simultaneous fit of as{Mzo) and ccq using resummed 0(Q;g)+NLLA QCD predictions with 
the ln(R)-matching combined with power corrections. The dotted lines represent an 
extrapolation of the fit result. 
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Figure 7: The energy dependence of (1 — T), (Af^), (Bt), (i?w), (C) and {y^) is shown. 
The sohd curves are the results of fits using perturbative 0{al) QCD calculations com- 
bined with power corrections while the dashed lines indicate the contribution from the 
perturbative prediction only. 
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Figure 8: Results for as(Mzo) and ao(2 GeV) from fits perturbative QCD predictions 
combined with power corrections to distributions a) or mean values b) of the event shape 
observables 1 — T, Mh or M^, Bic and i?w and C are shown. The error ellipses correspond 
to one standard deviation of the total error (38% CL) and take correlations from the fit 
and from systematic uncertainties into account as explained in the text. 
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